Self- Trapping of Bosons and Fermions in Optical Lattices 
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We theoretically investigate the enhanced localization of bosonic atoms by fermionic atoms in three- 
dimensional optical lattices and find a self-trapping of the bosons for attractive boson-fermion interaction. 
Because of this mutual interaction, the fermion orbitals are substantially squeezed, which results in a strong 
deformation of the effective potential for bosons. This effect is enhanced by an increasing bosonic filling factor 
leading to a large shift of the transition between the superfluid and the Mott-insulator phase. We find a nonlinear 
dependency of the critical potential depth on the boson-fermion interaction strength. The results, in general, 
demonstrate the important role of higher Bloch bands for the physics of attractively interacting quantum gas 
mixtures in optical lattices and are of direct relevance to recent experiments with 87 Rb - 40 K mixtures, where a 
large shift of the critical point has been found. 

PACS numbers: 03.75.Mn. 03.75.Lm, 03.75.Hh 



Ultracold multicomponent gases offer new insight into the 
atomic interspecies interaction. In optical lattices, the inter- 
play between interaction and tunneling is reflected in the com- 
plex phase diagram of boson-fermion mixtures. The multitude 
of predicted phases have been studied theoretically using var- 
ious approaches lUiSiiSSiiEiIl-In particular, 
the phase separation between bosons and fermions, the su- 
persolid phase |2l|3j], and the pairing of bosons and fermions 
forming phases of composite particles Q4Q have been investi- 
gated. Mixtures of bosonic 87 Rb and fermionic 40 K atoms in 
optical lattices were recently realized lfl2[ [l3ll and have drawn 
much attention due to an unexpected large shift of the bosonic 
phase transition between the superfluid phase and the Mott in- 
sulator lfl4l [l5ll . The effect, which is substantial even for a 
small ratio of fermionic to bosonic atoms, was controversially 
discussed Il3l[l6i[l7i[l8ll . and a relatively small influence of 
the boson-fermion interaction on the transition was proposed. 
In addition, the loss of coherence due to an adiabatic heating 
in the optical lattice was addressed B17lll8ll . 



In preference for a single-band Hubbard-type Hamilto- 
nian, the influence of interaction induced orbital changes were 
widely neglected in calculations performed for optical lattices. 
In this Letter, we show by means of exact diagonalization that 
the attractive interaction between bosons and fermions causes 
substantially modified single-site densities. The nonlinearity 
of the interaction leads to a mutual trapping of 87 Rb and 40 K 
atoms in the centers of the wells. For high bosonic filling 
factors the fermion orbitals, i.e., the local atomic wave func- 
tions, are strongly squeezed, which leads to a considerable 
deformation of the effective potential for bosons. The latter 
causes a large shift of the bosonic Mott transition towards 
shallower potentials which coincides with experimental re- 
sults Jl2, 13]. In this Letter, we show that the contribution of 
higher Bloch bands, usually neglected, can be very important 



for the physics of ultracold quantum gases in optical lattices. 

The interaction between fully spin-polarized neutral atoms 
in the ultracold regime can be described by a contact potential 
gS(r — r'). The strength of the repulsive interaction between 



87 



Rb atoms is given by g% = -^-ae, where as lOOao is 
the bosonic scattering length and ao the Bohr radius. The in- 
teraction between 87 Rb and 40 K atoms far from a Feshbach 
resonance is attractive with the strength (?bf = 
where n = " lBmp is the reduced mass and arr « — 205a 
B19I1 . The Hamiltonian of a Bose-Fermi mixture in an optical 
lattice is given by H = Hb + Hf + Hbf, where Hg describes 
the system of interacting bosons, H F the system of noninter- 
acting fermions, and H BP the interaction between bosons and 
fermions |Q1 01 ■ Using the bosonic and fermionic field opera- 
tors V>b(i") an d V"F( r )i me three parts of the Hamiltonian can 
be written as 
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with the periodic potential V(r) of the optical lattice. 

Numerically, we perform an exact diagonalization of the 
Hamiltonian ([TJ in a many-particle basis 12011 that includes 
higher orbital states and is truncated at a sufficiently high en- 
ergy (see lull for details). In our calculation, the bosonic and 
the fermionic subspace are diagonalized separately using a 
self-consistent interaction potential, which converges within 
a few cycles. For a known fermionic density pp(r) = (^p^f) 
the effective Hamiltonian of the bosonic subsystem is given by 
Hb S (pf) = Hb+9bf Jd 3 r p F (r) t/J|(r)^(r). The latter term 
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represents the interaction with the fermionic density and leads 
to a bosonic effective potential V^ n {pf) — V(r) + </bfPf(^)- 
Starting with the density of noninteracting fermions, we de- 
termine the boson density Pb{f) = (^Jj^b) by diagonal- 
ization of H^ 1 . Afterwards, a new fermion density can be 
calculated using the fermionic effective potential Vp eff (/3B) = 
V(r) + .9bf/0b(i 1 ). 

We discuss the experimental situation, where one fermion 
(tif = 1) and several bosons (1 < tib < 10) are present at 
each lattice site. To study the effect caused by the mutual in- 
teraction, it is instructive to restrict the particles to a single 
lattice site. Hence, we perform the diagonalization for atoms 
in a symmetric well with the shape Vb[sin 2 (fca;) + sin 2 (fcy) + 
sin 2 (fez)] [Eitl . To illustrate the orbital changes, we fit Gaus- 
sians to the calculated densities for bosons and fermions 
/ 27rcr B /F)~ 3 exp(' 
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PB/F{ r ) — "-b/f(v ^»^b/f; c^p^— / i^u^j-p, 
width of the Gaussians <t b /f a s a function of the bosonic fill- 
ing n B is shown in Fig. Q] for the lattice depths Vb = 8E R , 
14i? R , and 20E R , where E R — is the recoil energy of 
8 'Rb atoms. Since 40 K atoms are much lighter than 87 Rb 
atoms, the resulting density profile for a 40 K atom on a sin- 
gle lattice site is much broader than for the 87 Rb atoms. For a 
pure bosonic system (dashed lines), we observe an increase of 
the bosonic width with an increasing number of 87 Rb atoms, 
due to the boson-boson repulsion. The interaction with the 
single fermion leads to a compression of the boson density 
(solid lines), which becomes even narrower than the single- 
particle boson density. In fact, the bosonic width decreases 
slightly with an increasing and reaches a minimum at 
riB = 7 — 8. This effective attractive behavior is surprising, 
since the attractive interaction between bosons and fermions 
scales linearly with n B , whereas the repulsion of the bosons 
is proportional to n§. This effect is caused by the strong 
squeezing of the fermion density due to its effective poten- 
tial Vp B = V(r) + gspps(r), that is deepened linearly with 
the number of bosons (inset of Fig. |2). The increasing 
curvature of the fermionic effective potential, which equals 

J^y F eff | r=0 « 2V Q k 2 + V2^ 3 |ffB F | nn/oi using the Gaus- 
sian approach, causes the width of the fermion density to be 
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FIG. 1: (Color online) The width of the density profile for one 
fermionic 40 K atom (left) and ns bosonic 8r Rb atoms (right) for 
V = 8E R (+), V = UE R (o), and V = 20E R (x). The dashed 
lines are obtained for vanishing boson-fermion interaction. Because 
of the mutual attraction, the densities are substantially squeezed. The 
lattice constant a — n/k is 515 nm. 



similar to the bosonic one for n B = 3 — 4 and even narrower 
for riB > 4. In our calculation, the occupation of higher 
single-particle fermion orbitals, due to the squeezing of the 
fermion density, is roughly 10% for n B = 4 and reaches 40% 
for tib = 10. 

The effective potential V^ ff (pF) experienced by the 87 Rb 
atoms is plotted in Fig. [2] for Vq = 14E R . The curvature of 
the the bosonic effective potential ^ 5 -V [ | ff | r= o ~ 2Vb&; 2 + 

v2^ | <7bf I / <7p is strongly enhanced by the nonlinear depen- 
dence on the fermion width op, which leads to a deepening 
of the bosonic effective potential with increasing n B . The re- 
sulting energy gain in (Hbf) overcompensates the repulsion 
energy of the bosons. Instead of a stronger repulsion, an in- 
creasing 7Jb causes a self-trapping of the bosons in the deep- 
ened center of the effective potential. This effect is mediated 
by the squeezing of the fermion orbital due to the higher bo- 
son density. The compression of boson and fermion densities 
should be observable in the experiments by imaging a broad- 
ened momentum distribution or by using atomic clock shifts 
to measure the peak density directly J2J]. 

We now address the consequences of more than a single 
site but only one fermionic "impurity." We use the same di- 
agonalization method as above that includes the deforma- 
tion of orbitals for small quasi-one-dimensional chains rf2lll 
with 5 — 7 sites, a bosonic filling factor 1 < n B < 3, and 
15.Er < Vb < 20-Er. The self-trapping causes a local deepen- 
ing of the bosonic effective potential for the fermion occupied 
site. The gain in energy, due to the boson-fermion interaction, 
leads to a binding of several bosons to one fermion, which 
is a phenomenon resembling polaron physics 1 25 , 26A . Re- 
markably, we find that the fermionic impurity causes, for the 
above parameters, the localization of six additional bosons at 
its site. In experimental setups an additional confining poten- 
tial is used 112lll3ll . which establishes a finite atomic cloud, so 
that bosons and fermions occupy only the central lattice sites. 




FIG. 2: The effective potential V£ ff for nB bosonic 87 Rb atoms on 
one symmetric site created by one fermionic 40 K atom. The dashed 
line represents the unperturbed lattice potential V{x) for vanishing 
boson-fermion interaction. The inset shows the effective potential 
Vj? ff experienced by the fermion. 



3 



Because of the self-trapping, the bosons are attracted to sites 
that are occupied by fermions. If the fermionic filling np < 1 
or the fermionic cloud is smaller than the bosonic one (as in 
Ref. 11311 ). this leads to an increase of the bosonic filling fac- 
tor in mixtures in comparison with pure bosonic systems. In 
experiments, the predicted high bosonic site occupation will 
be accompanied by high losses due to three-body recombina- 
tions. 

In the following, we discuss commensurately filled macro- 
scopic cubic lattices and the implications on the Mott insu- 
lator phase transition that result from the deformation of the 
effective potential. In optical lattices, the bosonic superfluid 
to Mott insulator transition is triggered by the boson-boson 
interaction strength relative to the hopping amplitude. In the 
Bose-Hubbard framework 11141 Il5l 12711 . this is described by 
the ratio between the on-site interaction U and the hopping pa- 
rameter J that both depend on the depth of the lattice potential 
Vq. For boson-fermion mixtures, both U and J also depend 
explicitly on the filling factor. However, to estimate how the 
the interaction with the fermions influences the phase transi- 
tion, we use an effective Bose-Hubbard model, where the lat- 
tice potential is substituted by the bosonic effective potential 
V^ {t (ris) for a specific static filling factor tib- This approach 
leads to a renormalization of the parameters U and J in the ef- 
fective bosonic system. The on-site interaction U is obtained 
directly from V£ s (tib) and the hopping J by band-structure 
calculations using a finite periodic continuation of V^ (tib) 
1 28I1 . Qualitatively, the bosonic effective potential, which is 
shown in Fig. |2j reveals two important aspects: First, even 
for a single boson the minimum is deepened substantially and 
decreases further with an increasing number of 87 Rb atoms; 
second, the shape of the effective potential deviates strongly 
from the Vq sin 2 (kx) potential, leading to a broader barrier be- 
tween neighboring lattice sites. Both effects lead to a reduced 
hopping between neighboring sites, whereas the on-site inter- 
action energy is mainly increased by the stronger curvature of 
the effective potential. The decrease in J, which is larger than 
the effect on U (inset of Fig. [3J, causes narrower bands. 

The ratio U/J of the effective potential VJf ff (7iB) is plot- 
ted in Fig. [3] as a function of the lattice depth Vq, where 
the dashed black line corresponds to a pure bosonic system. 
The interacting system is represented by solid lines which are 
split for different filling factors ub, since the deformation of 
the effective potential grows with an increasing filling fac- 
tor riB- Even for = 1 the self-trapping causes a large 
shift of U/J (Vq) towards higher values, which is further en- 
hanced with an increasing ub- In the following, we use the 
renormalized values of U /J and known results for the criti- 
cal point (U/J) c to calculate the shift of the critical poten- 
tial depth V C . Following the mean field results, the critical 
ratio depends on the filling factor and obeys the relation 
(U I J) c = z\2nB + l + 2y/riB(nB + 1)], where z = 6 for a cu- 
bic lattice Q290 . In Fig. [3] the critical values for different filling 
factors ?iB are depicted as horizontal grey lines. The critical 
lattice depths V c (nB) for the pure bosonic system are given 
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FIG. 3: (Color online) Calculated ratio U/J plotted against the lat- 
tice depth Vo for filling factors he — 1 to tib = 10 and np = 1 (solid 
lines). The dashed black line represents a pure bosonic system. The 
boson-fermion interaction causes a large shift of the critical potential 
depth Vo - The inset shows U and J separately for the filling factors 
riB = 1, 5, 10 and for vanishing interaction. 



by the intersections of the dashed black line with the horizon- 
tal lines (open circles). According to the relation above, the 
phase transition for a pure bosonic system shifts to deeper lat- 
tices for an increasing filling factor n B . For boson-fermion 
mixtures, the intersections representing the critical potential 
depths are indicated by solid circles. In comparison with the 
pure bosonic system, the phase transition in mixtures is shifted 
substantially towards shallower lattices. 

The parameters in our calculations have been chosen such 
that a comparison with the experiment in Ref. II 1311 is possi- 
ble. The critical potential depth can be estimated from Fig. [3] 
as 12.3-Er for ub = 1, increases to 13.8i?R for rig = 4, and 
decreases slightly for higher fillings. Comparing the boson- 
fermion mixture and the pure bosonic system with the same 
bosonic filling, the shift AV C of the phase transition is given 
by the difference of the respective critical lattice depths (see 
arrow for ub = 2). The expected shift increases substantially 
with the filling factor, e.g., AV ° — —2.2Er for tib = 1 and 
AV C = -5.2£r for n B = 4. In Refs. EEl, pure bosonic 
systems are compared to Bose-Fermi mixtures with the same 
number of bosonic atoms, not accounting for the increase of 
the local bosonic filling factor due to the self-trapping, as dis- 
cussed above. In addition, the confinement in experiments 
causes the formation of shells with different filling factors n^. 
Therefore, the reported shift in Ref. lfl3ll of approximately 
— 5-Er is an average value of a system, where at the center of 
the lattice > 5 and np = 1. Accounting for the inhomoge- 
neous filling, we average our results for ub < 5 to ub < 7 
leading to a shift between — 4.2Er and — 5.1-Er, which is 
in good agreement with the experiment. Using atomic clock 
shifts B24I1 . the local filling factors could be determined, which 
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FIG. 4: The shift of the superfiuid to Mott insulator phase transition 
AVq in dependence on the scattering length ratio cief/ob for the 
filling factors tib = 1 to tie — 5. The dashed lines are obtained by 
keeping the fermion orbital fixed. 



tures in optical lattices. We found a self-trapping behavior of 
the bosons in their effective potential mediated by the inter- 
action with the fermions. Using a model with effective Bose- 
Hubbard parameters U and J, the expected shift of the critical 
potential depth separating the superfiuid phase from the Mott 
insulator was estimated. Our results reveal a strong depen- 
dence of the shift on the boson-fermion interaction strength 
and the bosonic filling factor, which could explain a large shift 
as reported in Refs. B12LI13I1 . and are applicable to future ex- 
periments using Feshbach resonances to tune the interspecies 
interaction. Theoretically, the full inclusion of orbital changes 
is a challenge for the efficient calculation of lattice systems. 

We thank F. Deuretzbacher for helpful discussions and K. 
Patton for reading the manuscript. K. B. thanks EPSRC for 
financial support in Grant No. EP/E036473/1. 



would allow a more accurate comparison of experiment and 
theory. 

Uniquely, experiments with ultracold atoms allow the pre- 
cise tuning of the interaction strength by Feshbach resonances. 
Hence, the shift of the critical potential depth AVq can also 
be studied in dependence on the scattering length obf between 
bosons and fermions, while the bosonic scattering length 
<2b = lOOao is kept constant. Such an experiment would al- 
low a closer investigation of the mutual interaction and the 
effects due to orbital changes. In Fig. [4] we present the calcu- 
lated shifts AUq for bosonic filling factors n B = 1 to n B = 5, 
where the solid lines are calculated allowing for orbital de- 
formations and the dashed lines by using a rigid fermion or- 
bital obtained by Vp eff (r) = V(r). The latter corresponds to 
a single-band approach for fermions, leaving the orbital de- 
grees of freedom only to the bosonic subsystem. Because of 
a much weaker deformation of the bosonic effective potential, 
both U and J are less affected than in the previous discussion 
i30ll . The scattering length obf, which enters linearly in the 
Hamiltonian (HJ, leads to an almost linear shift AV C . In great 
contrast, the dependence on <jbf for the self-consistent calcu- 
lation (solid lines), which fully includes the orbital degrees 
of freedom, is superlinear due to the self-trapping. Thus, the 
mutual deformation of the effective potential is enhanced by 
a larger scattering length as well as by a higher bosonic fill- 
ing factor as discussed above. Therefore, the assumption of 
a fixed fermion orbital is only suited for weak interaction and 
low bosonic filling. 

For repulsive interaction between bosons and fermions 
(obf > 0) the shift AVq becomes positive but remains rather 
small with a sublinear dependence on obf- In this region, how- 
ever, the separation of bosons and fermions plays a role, which 
is not accounted for in our calculations. 

In conclusion, we have shown that orbital changes in at- 
tractive Bose-Fermi mixtures ( 87 Rb - 40 K) are nonnegligible 
as they lead to a substantial deformation of the effective poten- 
tial and a squeezing of the effective orbitals. The results are 
therefore of fundamental importance for quantum gas mix- 
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